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We use a coherent Bragg diffraction method to impart an external momentum to ultracold bosonic 
atoms trapped in a one-dimensional optical lattice. This method is based on the application ol a 
single light pulse, with conditions where scattering of photons can be resonantly amplified by the 
atomic density grating. An oscillatory behavior of the momentum distribution resulting from the 
time evolution in the lattice potential is then observed. By measuring the oscillating frequencies, 
we extract multi-band energy structures of single-particle excitations with zero pseudo-momentum 
transfer for a wide range of lattice depths. The excitation energy structures reveal the interaction 
effect through the whole range of lattice depth. 

PACS numbers: 32.80.Rm, 34.20.Cf, 42. 25. Fx 



I. INTRODUCTION 

Cold atomic gases loaded in optical lattices have drawn 
great attention in the past few years as a versatile and 
powerful experimental system to study strongly corre- 
lated many-body problems By taking advantage of 
great controllability and perfectness, this novel type of 
matter features a promising candidate to achieve some 
important modeling systems and help tackling many un- 
addressed questions in multi-disciplinary fields of physics. 
As an example, important progress has been made in 
manipulation and characterization of interacting bosons 
loaded in an optical lattice, both for the low-lattice/low- 
temperature Bose-Einstein condensate (BEC) regime 
and for the high-lattice/high-temperature Mott Insula- 
tor (MI) or normal regime 0t3- 

Among these efforts, the understanding of excitations 
gives crucial information about the underlying system, 
for its ability to characterize the system's response in dif- 
ferent regimes. For shallow lattices, the gas behaves as 
superfluid and the excitation spectrum can be described 
by the Bogoliubov theory within a mean-field description 
of interaction effect. For deep lattices, the gas enters the 
strongly correlated MI phase, where the excitation ex- 
hibits a gap at low energies. Experimental characteri- 
zation of these features include the measurement of first 
band effective mass by Bloch oscillation the demon- 
stration of gap in the MI regime Q , and the momentum- 
resolved probe of multi-band dispersion by Bragg spec- 
troscopy [iMl • 

In this manuscript, we report a novel experimental 
scheme to detect multi-band excitation with zero quasi- 
momentum transfer for interacting bosons loaded in a 
one-dimensional (ID) optical lattice across a wide range 
of lattice depth. The elementary excitations are created 
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FIG. 1; (Color online) Schematic setup (top) and time se- 
quence (bottom) of our experimental scheme. After loading 
the BEC in a ID optical lattice with an exponential ramp 
(40 ms), a laser beam of 5/is is sent at 6^ = 24° from the lon- 
gitudinal axis, and is scattered into a symmetrical direction 
according to the Bragg scattering condition. After holding 
the lattice potential for another duration of r, both the lat- 
tice and the magnetic trap are abruptly switched off and an 
absorption image is taken after a time-of-flight (TOP) of 30 
ms. 



using Bragg scattering, where light waves are scattered 
collectively in accommodation to the periodic distribu- 
tion of bosonic atoms. As a response to this coherent 
process, the atoms will recoil with a momentum in unit 
of the reciprocal lattice, and can be excited to higher en- 
ergy bands with the same pseudo-momentum. By moni- 
toring the time evolution of the momentum distribution, 
we can extract information about the excitation energy. 
In shallow lattices where gases are superfluid, we can get 
the {n = 1) as well as the {n = 2) band excitation en- 
ergy. For deep lattices where the ground state is a MI, the 
large band gap prohibits dramatic population to higher 
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band, such that only the (n = 1) band excitation energy 
is obtained. Throughout the whole range of lattice depth 
variation, we also observe a sizable signature of interac- 
tion effect, due to the large occupation number for each 
lattice site. 

We stress that our scheme is different from Bragg spec- 
troscopy, where atoms are excited by a direct transfer of 
energy and momentum via a two-photon process stim- 
ulated by two distinct laser beams [l(]| - [l5| . In Bragg 
spectroscopy, the two laser beams form an effective one- 
dimensional optical lattice moving along a specified di- 
rection. If we see from the lattice frame, the atoms which 
are stationary in the lab frame are moving along the op- 
posite direction, and the Bragg spectroscopy scheme cor- 
responds to an effective diffraction of atomic wave func- 
tion on this optical grating. Therefore, this procedure 
usually requires a coherence condition for atoms and a 
beam duration of milliseconds to achieve distinct diffrac- 
tion peaks, which corresponds to significant population 
of specific excited states. 

In contrast, our experimental scheme involves only one 
incident beam, and the periodic density distribution of 
atoms in optical lattices are considered as gratings, from 
which the photons are scattered and cooperatively am- 
plified given a certain geometry to fulfill the Bragg scat- 
tering condition. Therefore, this method does not rely on 
phase coherence of atoms in different sites, and can be di- 
rectly applied to the Ml/normal regime. Besides, as the 
speed of photon is c, the building up time of scattering 
light field is negligible and this optically coherent pro- 
cess requires a much shorter light beam duration only in 
scale of microseconds, hence brings less heating or other 
complications to the atomic system. 

Our paper is organized as follows. In Sec. |ll] we de- 
scribe the experimental scheme based on the application 
of a single Bragg diffraction pulse on a gas loaded in a one 
dimensional lattice. In Sec. Illll we present our results on 
the oscillatory behavior of the momentum distribution 
observed as a function of keep time in the optical lattice 
after applying the Bragg pulse. Sec. IIVI is devoted to a 
more careful analysis of Bragg diffraction obtained with 
a single light pulse and used as a tool throughout this ex- 
periment. This diffraction is revisited with semi-classical 
theory and an interpretation in terms of collective light 
scattering; In Sec. |Vl we summarize our work and draw 
conclusion. 



II. EXPERIMENTAL DESCRIPTION 

The experimental setup and protocol is illustrated 
in Fig.l. We first prepare a BEC of about 2 x 10'^ 
*^Rb atoms in a QUIC trap, with longitudinal length 
L = 100 //m and transverse length I = 10 /^m. The 
BEC is then loaded into a ID optical lattice Vo\{x) = 
V\i -\- cos(2fcia;)]/2 along its axial direction by an adia- 
batic ramping of 40 ms. Here, V is the lattice depth, and 
Ul = 27r/AL is the reciprocal lattice spacing with Al the 
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FIG. 2: (Color online) (a) TOF images for V = 12Er. From 
top to bottom, the holding time r increases from zero with a 
step of 12/is. An oscillatory behavior of the first order peaks 
at ±2hkL can be clearly observed, (b) Similar results for 
V — 2)7 Er. Although the interference pattern is wiped out 
and no peaks is present, the oscillation of the momentum 
distribution can be still determined. 



lattice light wavelength. After holding the lattice for 50 
ms, we shine the Bragg pulse from a certain direction 9 
(see Fig. 1). The beam is red-detuned hy 5 = 1.3 GHz 
from the (5S'i/2,F = 2) ^ (5P3/2,i^' = 3) transition, 
and for a typical duration of At = t2 — ti — b fis. After 
the Bragg pulse, we further hold the system for a certain 
amount of delay time r, and suddenly shut down all the 
trapping potentials to get time-of-flight (TOF) images 
which are taken after 30 ms of free expansion process. 

Due to the presence of the optical lattice, the atoms 
are initially distributed into an array of quasi-two- 
dimensional slices, which are separated equidistantly by 
the lattice spacing Al/2 = 426 nm. For each slice of 
atoms, the Bragg pulse incident from the direction of 9 
will be reflected to the direction of {tt — 9), such that the 
optical lengths of different light paths are the same. Be- 
sides, the coherent condition for diffraction also requires 
the optical length difference between light paths scat- 
tered from adjacent slices is a multiple of the Bragg beam 
wavelength Xb- In our case of As = 780 nm, these co- 
herent conditions single out a specific incident direction 
of the Bragg pulse with cos 9 = Xb/^l, only at which 
significant collective scattering of light could occur. 

As a consequence of cooperative scattering (cf. Sec. 
lIVp . the photons would exert a recoil momentum oi2hkL 
to the atoms, leading to the occupation of different states 
with various bands but zero quasi momentum. After 
that, the momentum distribution experience an obvious 
oscillation due to the presence of optical lattices. This 
oscillation can be observed from a series TOF images 
by varying the holding time r. One typical result for 
V = 12i?r{ is shown in Fig. 2(a), where the system is 
in the BEC regime and a clear oscillation of the occupa- 
tion number of particles in p = ±2hkL is observed. In 
Fig. 2(b), we also present a similar result for V = ST-Br, 
where the high lattice potential destroys the phase co- 
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hcrcncc and the system becomes a normal gas. Notice 
that although the interference patterns no longer exist in 
the TOF images, the momentum oscillation can still be 
clearly determined. 



III. RESULTS AND ANALYSIS 

For shallow lattices where the system is in the BEC 
phase, atoms arc macroscopically distributed in the low- 
est quantum state which is described by the Gross- 
Pitaevskii (GP) equation 

2n2 



+ K>i(^) + + f/o|«'(^)|' 

2m 



(1) 

where T4xt(a;) = muj'^x^ /2 is the external harmonic trap, 
C/q = A-Kh^as/m is the interaction with as the s-wave 
scattering length, and ^ is the chemical potential. In the 
mean field level, the repulsive interaction acts as an effec- 
tive potential C/o'^o(2^) with tiq the condensate fraction, 
which tends to reduce the height of the optical lattice. 
This combined effective lattice Vcf[[x) = Vo\{x) + UonQ{x) 
has the same period as the original one, and forms a band 
structure for single particle state. The presence of BEC 
ensures that only the state at the band bottom is macro- 
scopically occupied. This state is characterized by the 
band index [n ~ 0) and pseudo- momentum ((7 = 0), de- 
noted by 



^{x,ti) = 4>n=o,q=o{x) 



R 



Wn^X^R), (2) 



where (pn.qix) is the Bloch state, Wn{r) is the Wannier 
function, and the summation over R runs over all lat- 
tice sites. After the Bragg pulse, an external momentum 
of 2hk]^ is transferred to the atoms. As this momen- 
tum transfer is resonant with the reciprocal lattice, there 
should be no new pseudo-mementum state being excited, 
and the system is described by 

'^{X,t2) = y^^an4'n,q=o{x) 

n 

= L^'/'Y.'^r.ci.^oe'''''^^, (3) 

n,£ 

where q = {2lhkL\(f)n,q) is the projection of the Bloch 
state to the plane wave state with momentum 2ihkL, 
and a„ is determined by the intensity and duration of 
the Bragg beam. Here, we should emphasize that the 
Bloch states 4>n,q{x) with g = are periodic functions of 
Ai/2, hence can be further expanded by a series of plane 
waves with t being integers. 

After the Bragg pulse, the wavefunction ^"(2;, t) evolves 
in the presence of the optical lattice as 



i2£hkL X 



Y,ayy^-\ (4) 



Here, we have dropped the subscript (g = 0) for simplify 
notation, and £"„ is the corresponding eigenenergy of the 
n"^ band. In realistic conditions, the summation over 
n is restricted to 71 < 2 due to negligible population of 
higher bands. Therefore, the average momentum (p) = 
(4'(t)|p|\E'(t)) and the atom number in the momentum 
mode Ni would acquire oscillatory behavior 



2hkL 



Cio cos(a;ioT + ipio) 

+ C21C0S{U!21T + ip2l), (5) 

D{q cos(i:jiot 4- ^10) + D^20 cos(w2ot + ^20) 

+D21C0S{UJ21T + (p2l), (6) 



where the three frequencies correspond to the energy 
gaps of cjio = (El - Eo)/h, W20 = {E2 ~ Eo)/h, and 
W21 = {E2 — Ei)/h, respectively, as illustrated in Fig. 
3(c). Notice that the frequency UJ20 = {E2 — EQ)/h is 
not present in the oscillation of (p). This is because the 
momentum distribution of states in bands (n = 0) and 
(71 = 2) are both even symmetric, hence make no net 
contribution to {p). The coefficients C's and D's are de- 
termined by 



2|a„a„| y ^ €c„c„. 



2|a, 



(7) 

the initial phase factors are determined by ipnm = 
angle(a^a„), the TVtot is the total atom number and the 
time average of Nn is = Y,n knPlcfJ^iVtot- 

Given the expressions of Eqs ([5]) and ([6]), we can ex- 
tract the three frequencies from evolution of the average 
momentum (p) [Fig. 3(a)], and the relative number of 
particles in p = ~2hkL state [Fig. 3(b)] via a fitting 
process, as depicted with solid (red) curves. One should 
notice that since the frequency W21 is relatively small, it 
would require longer evolution time r to be precisely de- 
termined. Considering the typical range of t < 250/iS in 
our experimental scheme, a more practical strategy is to 
fit the other two frequencies wio and 0^20 first, and rely 
on the relation u!2i — W20 — wio to get U!2i- 

In Fig. 4, we show the results of oscillating frequencies 
of LUiQ and W20 for various optical lattice depths ranging 
from 8 to 37 Ej^. Notice that the frequency 0^20 can only 
be obtained for lattice depth V < IdE^, as the strong 
lattice potential prohibits significant excitation to the 
(n = 2) band. For shallow lattices with V < 25Eii, the 
system is in the BEC regime presenting clear interference 
pattern in the TOF images, and the discussion above re- 
mains valid. If we neglect the interaction effect between 
atoms, the band structure can be calculated based on 
the bare optical lattice potential, and the mutual gaps 
between different bands are indicated as solid lines. As 
one would expect, this estimation is systematically above 
our experimental observation, since the repulsive interac- 
tion tends to reduce the lattice depth to a lower effective 
value. 
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FIG. 3: (Color online) From the TOF images, we can ex- 
tract (a) the average momentum of the cloud and (b) the 
relative occupation number of the —1 order with p — —2hkL- 
By fitting the results with Eqs. (|6]) and (O (red solid), we 
can extract the frequencies between multi energy bands, as 
schematically drawn in (c). The results correspond to the 
case oiV = 12Er, and all date points are taken as average of 
more than three experimental trails. 



One theory to incorporate the interaction effect is dis- 
cussed in Ref. [3, where Choi and Niu analyzes the 
interaction effect for the bottom (n = 0) band and 
gives an approximate expression for the effective poten- 
tial Vcs = 1^/(1-1- AC) by assuming a nearly uniform con- 
densate distribution. Here, C = nnoas/k'j^ is the dimen- 
sionless interaction strength and uq the average three- 
dimensional condensate density. By using the effective 
potential, the oscillating frequencies can be extracted as 
mutual gaps between different bands, as illustrated in 
Fig. 4 (dashed lines). Another method to capture the in- 
teraction effect is to use a tight-binding model , which 
is more reasonable for deeper lattices since the Wannier 
functions for all bands are well localized within one lat- 
tice site. Thus, we can calculate the bottom and first 
exited band structures via a variational method by ne- 
glecting higher order overlaps between sites, and extract 
the frequency wio (dotted Hne in Fig. 4). We emphasize 
that neither of the theories above has included finite tem- 
perature nor inhomogeneity effect, and their comparison 
to our experimental data can only be understood to a 
qualitative level. 
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FIG. 4: (Color online) Oscillating frequencies ujio (solid 
dots) and lu20 (hollow dots) in the unit of recoiled frequency 
i^R = Er/Il. In the BEC regime with V < 25En (left), ex- 
perimental results are systematically below the bare lattice 
calculation (solid), where interaction effect is not taken into 
account. To incorporate the interaction effect, an effective 
potential theory (dashed) and a tight-binding model (dotted) 
are considered. In the Ml/normal regime with V > 25-E'r 
(right), the oscillating frequency uio is compared with the 
trapping frequency of the lattice potential (dashed-dotted), 
where we use Us — 4.7 nm, and no — 10^^ cm~^. 



Up to now, we have restricted our discussion to the 
shallow lattices regime where the system is in the BEC 
phase. By further increasing the lattice height over 25 
En, the strong local confinement would eventually break 
the global phase coherence between different sites, and 
the system leave the BEC regime. In this case, since 
the recoil energy from diffracted photons Eb ~ 4:Efi is 
much smaller than the confining potential, atoms that 
arc kicked by the Bragg pulse are mainly localized within 
one lattice site and swing to and fro around the site cen- 
ter. If we neglect the interaction effect, the oscillating 
frequency is thus the trapping frequency of the lattice 
potential, as shown in Fig. 4 (dashed-dotted line). No- 
tice that this bare lattice estimation is significantly above 
the experimental measurements. A detailed description 
of the interaction effect would require a solution of the 
hydrodynamic equation, which is beyond the scope of the 
present manuscript. 



IV. SEMICLASSICAL THEORY FOR THE 
COOPERATIVE SCATTERING PROCESS 

Finally, we notice that the oscillating frequencies be- 
tween different energy bands can be identified only for 
fairly deep optical lattices with V > 8Ey{ . For shallower 
lattices, we could not observe significant oscillation. This 
can be easily understood by reminding that the incident 
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FIG. 5: (Color online) Average momentum of the cloud from 
TOF images upon released from the trap right after the Bragg 
pulse of 5 /iS. Each data point is an average of more than 
three experimental trials. This result can be well described 
by a semiclassical theory (solid, see text) , with fitting coupling 
constant g = 0.9 MHz. 



photons are coherently scattered and amplified by the 
atomic grating, which becomes more apparent in deeper 
lattice potentials. This scattering processes can be de- 
scribed by using a semiclassical theory [l^. We start 
with the MaxwcU-Schrodingcr equation, which describes 
the evolution of atomic wave function vE'(x, t) and the pos- 
itive and negative frequency components of the classical 
electric field E(±): 



(d ■ E(~h(d ■ E(+') 



= cV^E^i 



h6 



p(±). 



(9) 



the lattice light, Xs is the coordinate along the direction 
of the scattered beam. The second term of Equation 
(fTni 's right side represents that an atom at absorb 
a photo from the pumping Bragg beam and emit another 
into the scattered beam. During this process it receives 
recoiled momentum of 2hkL and is transmitted into 
mode. On the other hand, the third term shows an atom 
from "ff+i mode absorb a photo from the scattered beam 
and emit one into the pumping beam and is transferred 
backward into VPf mode. Equation PT|) shows the evolu- 
tion and propagation of the scattered beam. Note that it 
is the atomic grating J^e ^^'^e+i ^^^^ initially transfers 
photons into the scattered beam from the pumping light. 

Assuming that the amplitude of light field varies 
slowly, and neglecting the spatial distribution of ^'^ and 



[Tgl . we get: 
dc± 



- = ±3(c^ziCo + CoC±i)co, 



dt 

^ = .9(|5.ip-|cip)co, 



(12) 
(13) 



where Ci{t) = {2ihkL\^{t)) , the amplitude of mo- 
mentum mode, can be solved numerically from the above 
equations, together with the initial conditions 5^(0) = 
<4i=o- 9 = \di2\'^\£p\'^i^snL/ {2£och^S'^) represent the cou- 
pling of the light with BEC of length L. Using the rela- 
tion ci = ^„ c^fln we can determine the coefficients an 
and their evolution, hence obtain the oscillation process 
discussed above. In Fig. 5, we show the average mo- 
mentum {p) ~ 2hkL J2e^ ' l^^l^ cloud right after a 
Bragg pulse of 5/iS duration, and compare with the semi- 
classical prediction. For shallow lattices, the momentum 
change is relatively small such that the oscillatory behav- 
ior upon time evolution is blurred. However, for deeper 
lattice potentials, significant amount of momentum can 
be transferred to atoms via this coherent scattering pro- 
cess. 



where H ~ ~5I7^^ + ^^i^) the single atom Hamilto- 
nian, S the detuning of the laser from the electronic 
transition, d the atomic dipole moment, and F^^^ = 
— ' ' ' 'E^^' the polarization. By decomposing the 
atomic wave and electric field as 4* = J^e ^'^e*'^^'^^^, and 

g(±) _ g±g-i(wp-fcia;-fej2!) 

get: 



^^^ — i{iAJs-\-kLX — kzz) 



d_ 

dt 



(14 



+ |2 



del 
dt 



+ c 



h^5 

det _ -ii|d|V 
dxs 2eoM 



(10) 



(El 



(11) 



where Ep and 



arc the amplitude of the pumping light 
and the scattered light, ujp and uJs are the correspond- 
ing frequencies, is the wave vector perpendicular to 



V. CONCLUSION 

In summary, an experimental investigation of multi- 
band energy spectrum for a bosonic system loaded in a 
ID optical lattice is demonstrated. A coherent Bragg 
light scattering process to excite atoms to higher energy 
bands with zero pseudo-momentum transfer, and oscil- 
latory behavior of average momentum upon holding the 
sample of atoms in the optical lattice for various evolu- 
tion times, arc first reported. We extract the oscillat- 
ing frequencies, which correspond to the energy gap be- 
tween different bands, for a wide range of lattice potential 
depth. Throughout the whole range of lattice depth, we 
observe the effect of repulsive interaction between parti- 
cles, which tends to reduce the lattice height to a lower 
effective value. Our method is a coherent process where 
photons are scattered and amplified by the atomic den- 
sity grating, while the previously reported Bragg spec- 
troscopy requires a high coherence of atoms. 
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